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Despite its remarkable agreement with gravity tests in the solar system, general relativity is 
challenged at larger scales. Large scale tests performed in the solar system by the Pioneer 10/11 
probes have failed to confirm the expected gravity laws. We discuss here extensions of general 
relativity which preserve its geometrical nature while showing an ability to conciliate the Pioneer 
anomalies with other gravity tests. Further testable phenomenological consequences are also briefly 
discussed. 
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The regularly improving agreement of gravity tests performed in the solar system with general relativity (GR) [1, 2] 
has imposed the latter as the standard theory for gravitation. However, anomalies appear at very large scales, as in 
the rotation cuves of galaxies [3, 4] or in the relation betweeen redshifts and luminosities of some type II supernovae 
\& • [5, 6]. They can be cured either by introducing "dark matter" and "dark energy" or by modifying the gravity law at 
1 \ large scales [7-9]. 

Testing gravity laws at large scales is clearly of great interest in this context. Such a test was performed in the 
q ' extended mission of Pioneer 10/11 probes and it failed to confirm the known laws of gravity [10, 11]. The Doppler 
O^, tracking data of the Pioneer probes exhibited deviations from the predictions of GR, the Doppler velocity residual 
varying linearly with time. This "Pioneer anomaly" is thus naturally represented as an acceleration ap directed 
bi), towards the Sun, with an approximately constant amplitude over a large range of heliocentric distances (AU = 
l— ~~ '■ astronomical unit) 

a P = (0.87 ± 0.13) nm s~ 2 , 20 AU < r < 70 AU (1) 

m ; 

The anomaly adds to the puzzle of large scale gravity, as it may reveal an anomalous behaviour of gravity already at 
\Q . scales of the order of the size of the solar system. 

Although a number of different mechanisms has been considered to this aim [12-15], the anomaly has escaped up 
to now all attempts of explanation as a systematic effect generated by the spacecraft itself or its environment [16, 17]. 
This has motivated efforts devoted to recovering and reanalyzing Pioneer data over the whole duration of the Pioneer 
10/11 missions [18], as well as proposals for new space missions dedicated to the study of deep space gravity [19]. 

The present paper is devoted to the key issue of a confrontation of the Pioneer anomaly with theories of gravity. We 
will show that a family of metric extensions of GR can be defined in a natural way and that it shows the capability 
• »"H , of making the Pioneer anomaly compatible with other gravity tests. We will also briefly discuss further testable 
' phenomenological consequences of this framework. 
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II. GRAVITY TESTS IN THE SOLAR SYSTEM 



We first briefly review the agreement of GR with tests performed up to now in the solar system. Before discussing 
this question, one must first recall the two building blocks which found GR as a gravitation theory. The first one, 
the equivalence principle, states the universality of free fall and gives GR its geometric nature. Violations of the 
equivalence principle are usually looked for by following the difference r\ in the free fall accelerations undergone by 
two test bodies having different matter compositions. This parameter r\ is constrained by modern experiments to 
remain extremally small, below the 10 -12 level, so that the equivalence principle is one of the best tested properties 
of nature. The validity of the equivalence principle has also been tested over a large range of distances, from the 
millimeter [20] to the sizes of Earth-Moon [21] or Sun-Mars orbit [22, 23]. Obviously, the level of precision attained 
by tests precludes an explanation of the Pioneer anomalies as a violation of the equivalence principle, as those would 
occur at the 10 -3 level (at the largest explored distances). 

The second basic property characterizes GR as a field theory by prescribing a particular coupling between gravitation 
and its sources. This coupling is equivalent to the gravitation equations which allow one to determine the metric 
tensor from the distribution of energy and momentum in space-time. According to GR, the coupling between curvature 
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tensor and stress tensor takes a particularly simple form [24-26], as one curvature tensor (Einstein curvature E^ v ) is 
simply proportional to the stress tensor T^ v . This coupling then involves a single constant, the Newton gravitation 
constant Gn 

— Rfj,u ~ ~^9^vR = ~~~ (2) 

Such a relation may be seen as the simplest way of relating curvature and energy-momentum tensors. As a result 
of Bianchi identities, E^ u corresponds to the particular curvature tensor having a null covariant divergence, like the 
energy-momentum tensor. Whilst the first property follows from Riemannian geometry, the second one expresses the 
conservation of energy and momentum and describes the geodetic motion undergone by freely falling test bodies. One 
notes however that these properties do not by themselves imply Einstein-Hilbert choice. 

In standard gravitation theory, the metric tensor in the solar system is deduced by solving Einstein-Hilbert equation 
(2) with the Sun as gravitation source. Rotation and non sphericity of the Sun may be discarded for our present 
purpose, so that the gravity source is described as a point-like motionless body of mass M and the metric may simply 
be written in terms of the Newton potential (f>. The solution is conveniently written in terms of spherical coordinates 
(c denotes light velocity, t and r the time and radius, 9 and <p the colatitude and azimuth angles) with the gauge 
convention of isotropic spatial coordinates 



ds 2 = g oc 2 dt 2 + g rr (dr 2 + r 2 (d9 2 + sin 2 9d(p 2 )) 

+ 20 4 
G N M 



.goo = 1 + 20 + 20 2 + ... , 5rr = -l + 20 



H « 1 (3) 

err 

Then GR, or Einstein-Hilbert equations, (2) can be tested by comparing the gravitation effects entailed by solutions 
to this equation with those predicted by alternative theories. For this purpose, one embeds GR within a family 
of well defined and computable theories. This is usually done by confronting GR with the family of parametrized 
post-Newtonian (PPN) metric theories [27-31] 

.g oo = l + 2a0 + 2/30 2 + ... , g rr = -1 + 2 7 + . . . (4) 

Among the parameters that can be considered, the constants a, (3 and 7 appear as the most relevent ones. The first 
one can be set to unity by redefining Newton constant Gn, GR corresponding to 7 = (3 = 1. The parameters 7 or (3 
affect geodesic motions in a way which can be evaluated, leading to comparison with observations and providing best 
fit values for these parameters (4). 

Experiments performed up to now have put strict bounds on the potential deviations 7 — 1 and (3—1 which have 
improved with time. From Dopplcr ranging on Viking probes in the vicinity of Mars [22] to deflection measurements 
using VLBI astrometry [32] or radar ranging on the Cassini probe [33], possible values of 7 — 1 have decreased to 
less than a few 10~ 5 . Precession of planet pcrihelions [34] and polarization by the Sun of the Moon orbit around the 
Earth [35] constrain linear superpositions of (3 and 7 fixing (3 — 1 to be smaller than a few 10 -4 . 

A complementary test of GR is to check the dependence of the Newton potential 4> (component goo in (3)) on the 
radial coordinate. Modifications of Newton standard expression are usually parametrized in terms of an additional 
Yukawa potential depending on two parameters, the range A and the amplitude a measured with respect to Newton 
potential <p 

5<p(r) — ae~^4>(r) (5) 

A Yukawa correction has been looked for at various distances ranging from the millimeter [20] to the size of planetary 
orbits [36]. The accuracy of short range tests has been recently improved, as dedicated experiments were pushed to 
smaller distances [37], while accounting for Casimir forces which become dominant at such ranges [38]. On the other 
hand, long range tests of the Newton law, performed by monitoring the motions of planets or probes in the solar 
system, also agree with GR with good accuracy, for ranges of the order of the Earth-Moon [21] to Sun-Mars distances 
[22, 39, 40]. But these results still show [36, 41] that windows remain open for violations of the standard form of 
Newton force law at short ranges, below the millimeter, as well as long ones, of the order of or larger than the size of 
the solar system. 

Thus, although they confirm the metric character of gravitation and provide strong evidence in favor of a metric 
theory lying very close to GR, gravity tests still leave room for alternative metric theories, at ranges which have been 
well tested. They furthermore point to a potential scale dependence of violations of GR. Such extensions of GR may 
in fact arise naturally, in particular when effects of radiative corrections are taken into account [42-46] . 
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III. METRIC EXTENSIONS OF GR 



After the preceding discussion, we focus our attention on metric theories of gravitation while releasing the assump- 
tion that the solar system is exactly described by the metric calculated from GR. In other words, we consider metric 
extensions of GR, that is to say metric theories of gravitation lying in the vicinity of GR. 

Such extensions are more easily explained in their linearized approximation where the gravitation field is represented 
as a small perturbation h^ u of Minkowski metric r}^ 

V, = diag(l, -1,-1,-1) , | VI « 1 (6) 

In the linearized approach, the field may equivalcntly be written as a function of position x in spacetime or of 
a wavevector k in Fourier space, so that Riemann, Ricci, scalar and Einstein curvatures simply read in momentum 
representation and at first order in h^ v 



k\k y h^p k\kph^i/ — k^kyhxp + kpkph\ v 



R 
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Rfiu R fi\v 7 R R^ 1 i E^ v R^iv t]^v 2 (^) 

Then, gravitation equations which generalize Einstein-Hilbert equations (2) may be written under the form of a linear 
response relation between Einstein curvature and the energy-momentum tensor [44] 

E»M =X^x P [k] T x P[k] , fc^A P [fc]=0 (8) 

GR (2) corresponds to a simple local expression for the linear response function [24, 25]. Generalized linear response 
equations (8) may arise for example from radiative corrections of GR through the coupling of gravity to quantum fields 
[42-46]. Then, as gravitation couples to fields with different conformal weights, the resulting gravitation coupling 
differs in the corresponding sectors. The integral equation (8) takes a simpler form when the energy- momentum tensor 
of the source is static and pointlikc (M is a mass localized at the origin) 

= StrtSrfToo, Too = Mc 2 S(k ) (9) 

In their linearized approximation, extended gravitation equations (8) are then solved for the static point-like source 
(9), using projectors on the two sectors of different conformal weights 

E - E {a) + E {1) ir = n - 

^IlJ = g — I — — J oo, E^J = — loo 

G (0) = G N + SG (0) , G (1) = G N + <5G (1) (10) 



The generalized gravitation equations may then be written in terms of two gravitation constants which differ in the 
two sectors and depend on momentum or length scale. Two running coupling constants G^ and G^ thus replace 
the unique Newton gravitation constant G^, introducing a dependence on the conformal weight as well as on the 
length scale [47, 48]. 

Generally, as we restrict our attention to extensions lying in the vicinity of GR, deviations from GR can be 
considered as small anomalous terms and calculated perturbatively. As GR corresponds to Einstein curvatures 
vanishing everywhere except on gravity sources, that is the Sun in the solar system, this assumption allows for 
simplifications. For instance from Bianchi identities, transversality takes a simple form and Einstein tensor can 
still be decomposed on the sectors of traceless and traced tensors as in the linearized approximation [47, 48] (these 
expressions however preserve their full non linear dependence in Newton gravitation constant Gn) 

E=[E] GR + 6E, [E] GR = where T = 

SE = SE^ + SEW (11) 



4 



The information contained in the anomalous running coupling constants SG^ and 5G^ is equivalently encoded in 
the spatial variations of Einstein curvature SE^ and SE^ of a solution of gravitation equations, or also, of two 
independent components (for instance in spherical coordinates, 5E® and SE T r ). 

We begin with the solution of GR gravitation equations (2) for a static pointlikc source (9), which we write in 
isotropic coordinates (3) 



Metric extensions of GR may then be considered as perturbations of the standard solution (12) 

300 = [#Oo] GR + ^500, 9rr = [#rr] GR + Sg rr (13) 

The perturbed metric components (13) should then correspond to the anomalous Einstein curvatures (11) derived 
from the modified gravitation equations. Although not knowing the precise form of these equations, one may nonethe- 
less write a general form for their solutions. In particular, the two independent components of Einstein curvature 
characterizing the solution may be chosen so that they lead to simple expressions for the metric components 

8E^2u\8^ N -S^ P )'\ u = — - — ■ ()'=* [ 500 ] GR = l _ 2 -^ u 

SE; = 2u 3 6<S>' P (14) 

The metric lying in the vicinity of GR and having Einstein curvatures as given by equations (14) may then be expressed 
in terms of two anomalous potentials [49] 

<%>o _ „ f (S$n - 6<S> P y j , „ f 5<f> P ' ^ | [g 00 ] GR - 1 f (5$ N - 5$ P )' ^ 



= 2 fv*»- 2 **ry du+2 f. 



boo] GR J [ ff0 o] GR J M GR [ ffQ0 ]| R J [ goQ ^ R 

[9oo] GR J 3 du (15) 



j— — - -2 [goa\ GR / ; — -a 

,Srr| ™ J [5oo] GR 



At this stage, it is worth comparing generalizations obtained in this manner with the widely used PPN Ansatz [31]. 
This PPN Ansatz may be written as a particular case of the more general extension presented here 

S<S> N = ((3- l)^ 2 + 0(cf> 3 ), (5$ P = -( 7 - l)cf> + 0(<j> 2 ) 

5E° = ±0(4> 2 ), 6K = 4r(2( 7 - 1)0 + O(0 2 )) [PPN] (16) 

One remarks that the PPN Ansatz corresponds to non zero values for Einstein curvature bE T r apart from the source, 
at first order in GM , but not for 5Eq which vanishes at the same order. But, as seen later, the PPN Ansatz does 
not have the ability to account for the Pioneer anomaly, which is associated to a non vanishing curvature 6E®, with 
moreover a scale dependence which differs from that described by the PPN Ansatz. 



IV. PHENOMENOLOGICAL CONSEQUENCES 

We now discuss some phenomenological consequences of the metric extensions of GR. First, we show that they 
lead to effects on Doppler and ranging measurements which take the form of Pioneer like anomalies. To describe 
such observations, one must introduce the time delay function, defined as the time T taken by a light-like signal to 
propagate from a spatial position xi to another one X2 (representing the coordinates of two different points) 

T(x 1 ,x 2 ), x a = r o (sin0 a cos<p a ,sin0 a sinvj a ,cos0 o ), a =1,2 (17) 

The observable used in ranging corresponds to half the time elapsed on Earth between the emission time and the 
reception time of a signal exchanged between the station on Earth and the probe. Moreover, as the Pioneer probes 
were followed by Doppler velocimetry, their monitoring is more appropriately discussed in terms of the Doppler 
acceleration a, [11] which is just the second time derivative of the round trip time delay 

T(xr,x 2 ) + T(x+ x 2 ) a_d 2 At 

At = 2 ' - C = ^P ( 18 ) 
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In the case of a static isotropic space-time the time delay function only depends on three real variables, which can be 
chosen as the radial coordinates r\ and r 2 of the two points and the angle <f> between them as seen from the gravity 
source. Its exact form is obtained by solving Hamilton- Jacobi equation for a light ray [48, 49] 

(r)dr p pdr / r 2 



cT(r u r 2 ,d>)= g "° , cj>= / 



(r) — P— Jr-i / 9rr l r \ _ 2— 

goo y ' r 2 V 9°° *' 2 

cos0 = cos^i cos 6*2 + sin^i sin# 2 cos ((p 2 — <pi) (19) 

Perturbations of GR metric then entail perturbations of the time delay function 5T and of the probe trajectory 
(Earth motion xi is assumed to remain fixed) which are used for determining the round trip time delay (18). 
Both perturbations affect the observed Doppler acceleration (18), so that when comparing the latter with the result 
predicted by GR, the residual difference does not vanish and may be interpreted as an acceleration anomaly affecting 
the Doppler observable [50] 

^=P1gr + <""> x 2 = [x 2 ] GR + <5x 2 

a ee [a] GR + 5a (20) 

For an extension in the neighborhood of GR, the acceleration anomaly may be obtained in terms of the two metric 
components which characterize the perturbed static isotropic metric (13,15). The expression obtained for the accel- 
eration anomaly in the general case [50] simplifies in the practical case of the Pioneer probes, as several contributions 
may be neglected in a first step. First, the GR metric components differ from 1 for [ffoolcR an d from -1 for [sv-t-Jgr by 
the Newton potential, which at the distances explored by Pioneer probes is smaller than 10~ 9 . Corrections induced 
by the square of the ratio of the probe velocity to light velocity are of the same order (~ 10~ 9 ). Then, the impact 
parameter of the signal [/o] GR (~ ri) is at least 20 times smaller than the probe distance [r 2 ] GR . Finally, the probes 
follow almost linear trajectories. Neglecting all the corresponding correcting terms, the acceleration anomaly reduces 
to the sum of two main contributions [50] 

5a ~ 5a scc + 5a ann 



sec 

,2 



6a scc ~ -y<9 r (<5.9oo) + N GR j (5( - g0 ° grr - ) - <5 9oo j - C —d 2 r [g o] GR 5r 2 



The term <5a sec contains secular contributions proportional to metric anomalies and depending on the probe accelera- 
tion, as well as to a range ambiguity <5r 2 . As unfortunately no range capabilities were available for the Pioneer probes, 
the latter should be determined by fitting the probe trajectory with the Doppler signals using a given model. The 
term 5a ann is a modulated contribution, proportional to the relative angular velocity between the probe and Earth 
and to an anomaly of the impact parameter 5 p. Using a simplest form for the impact parameter [p] GR , modulated 
by Earth rotation, and taking as reference the acceleration on Earth due to the Sun gravitation field, the secular and 
modulated anomalies may be expressed in terms of the potential anomalies in the two different sectors [50] 

5a scc ~ -c 2 d r 5<S> N (r 2 ) + ai — \ + 5<5> P (r 2 ) + 25^ N (r 2 ) 

r 2 z {_ r 2 

dA 

5a ann ~ -Acos(2n 1 (t - iconj)) - sin (2^i(t - t conj )) 

r 2 [ r 2 J ri r z 

ai = ^lE = njn ~ 6 x 10- 3 ms- 2 (22) 
ri z 

The time t con j corresponds to the closest approach between Earth and probe. The anomaly <5a ann depends on an 
amplitude A which is determined by the range ambiguity and the two potential anomalies, and is modulated at twice 
the orbital frequency 

The secular and modulated anomalies involve different linear superpositions of the range ambiguity and the poten- 
tials (22). Although the modulated anomaly seems to be larger than the secular one by a factor r 2 /r\, this turns out 
not to be the case in practice [11]. Indeed, as the data analysis process is based on the a priori assumption that there 
is no modulated anomaly in the physical signal, the best trajectory fitting the data tends to produce a quasi null value 



6 



for the modulated anomaly, i.e. a nearly perfect compensation of its different contributions. Such a compensation has 
to be effective within a fraction of the order of ri/r 2 but cannot hold over a long period of time, due to the different 
range dependences of the compensating terms and also to the maneuvers which change the constants of motion, hence 
the range ambiguity. A precise estimation of the modulated anomaly would then require a complete solution of the 
motion including a detailed description of the maneuvers. The latter being unavailable, we shall use the previous 
description to obtain a rough quantitative estimate of the secular anomaly. Inserting for the range ambiguity in (22) 
the value minimizing the modulated anomaly and keeping leading contributions only, one obtains a general expression 
for the secular anomaly which only depends on the two potential anomalies 

5a sec ^-c 2 a r ^(r 2 ) + a 1 r 1 2 {A C 5Z P {r)% + 6 -^pl\ (23) 

L r 2 Jri r r 2 ) 

One should note on the final form of the anomalous acceleration (23) its non linear dependence on the two gravitation 
potentials, as it involves a contribution which is the product of the standard Newton acceleration (first sector) by 
the potential in the second sector. In particular, this result differs significantly from those obtained from preliminary 
studies using the same framework [48, 49] but neglecting the modulated anomalies, which in fact appear to play a 
crucial role. 

According to the expression (23) for the secular anomaly, a constant anomalous acceleration (5a scc — —ap/c 2 = 
— as observed in the distance ranges of the Pioneer 10/11 probe, may be obtained in several ways. First, assuming 
a perturbation limited to the sector of Newton potential, one obtains that the latter must depend linearly on the 
heliocentric distance 

5$n * j (24) 

However, the effects of a metric keeping the same form (24) between Earth and Mars could not have escaped detection 
in tests performed with martian probes [11, 39, 41]. The perturbation (24) would produce a change of planet orbital 
radii, leading to range variations of the order of 50km to ^lOOkm, while the Viking data constrain these measurements 
to agree with standard expectations at a level of ^100m to ~T50m. One should further note that the effect of the 
Shapiro time delay in the range evaluation has here a negligible influence [51]. A remaining possibility is that the 
perturbation (24) be cut off at heliocentric distances corresponding to inner planets [52], a possibility which should 
then be confronted with available data on the ephemerids of outer planets. 

If on the other hand the perturbation is limited to the second sector, a constant anomalous acceleration is obtained 
from a quadratic dependence of the potential SQp on the heliocentric distance 

„2 2 

"'"-mm (25) 

The metric perturbation (25) also has an influence on the already discussed range measurements on martian probes. 
If the dependence (25) were extended without modification to the inner part of the solar system, the Shapiro effect 
would lead to range variations which stand in slight conflict with Viking data. But this discrepancy may be cured 
by mildly cutting off the simple dependence (25) at the orbital radii of inner planets, without affecting predictions 
made for the Pioneer probes. In any way, a precise analysis of the anomalous acceleration (23) should consider the 
most general form depending on the potentials in both sectors (such as a linear combination of the two simple models 
(24) and (25)). The evaluation of correlated effects of anomalies in the two sectors will be necessary to confront 
observations with theory in a significant way. 

We now discuss effects of the extended metric in the inner part of the solar system. An important success of GR 
concerns its agreement with planet ephemerids. In the case of metric extensions of GR, and for the nearly circular 
orbits followed by planets, the anomalous precession may be given as an expansion in the orbit eccentricity e. The 
result, derived and discussed in [49], is just given here 

<5Atu „ c 2 u „ e 2 u 2 ( , , „ , -.a c 2 u 



2tt 



The well known expression of the anomalous perihelion precession in the PPN formalism [1] is recovered by inserting 
in (26) the particular expressions of 5&n and <5$p for the PPN metric (16). Keeping only the leading term, one 
obtains the anomalous perihelion precession in terms of anomalous parts of Eddington parameters (3 and 7 



SAzu\ PPN GM 



2tt 



(2( 7 -l) -(/3-1))— u (27) 
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The more general expression (26) may thus be thought of as a generalization of the anomalous PPN parameters (3—1 
and 7 — 1 to functions (<5<i>jv)" and (uSQp)" which may depend on the heliocentric distance. Expression (26) for 
perihelion precession anomalies may then be used to constrain the potential anomalies for ranges around a few AU 
where planet ephemerids are precisely known. 

Finally, we come to the discussion of light deflection experiments by the Sun gravitation field, which amount to 
explore the potentials $ n and <I>p in the Sun vicinity. To this aim, we first rewrite the definition of the anomalous 
part 5-y of the Eddington parameter. We decompose the relevant combination of the two potentials 2<5<I>Ar — (5$p and 
separate the part singular at r = from the regular part describing the long range effect 

2J$ JV (r)-<y*p(r) = -^^ + ^rCo(r) (28) 

The coefficient Go of the singular part in (28) corresponds to the standard PPN modification of Eddington parameter 
7. One specific feature characterizing the extended metric is then the dependence of the anomalous part of Eddington 
parameter ^7 on the impact parameter p (of the electromagnetic signal deflected by the Sun gravitation field), which 
is determined by the regular part of the anomalous potential 

*r(p) = -7T + ^ (29) 

This dependence, which allows one to reveal the long range variation of the potential 4>p, can be looked for by keeping 
trace of the whole functional dependence of the anomalous deflection angle during an occultation by the Sun 

G N M 8 ( Ar ir2 

dp y 5 ^— 



6Ae^-^-[6 7 (p)ln^) (30) 



For experiments within the solar system, r\ has to be replaced by the Sun-Earth distance and r 2 by the Sun-probe 
distance. When the emitter is a source far outside the solar system, r 2 has to be replaced by a cutoff distance. Evidence 
of a variation of 7 would provide a direct validation of the extended framework. This aim might be attainable through 
a reanalysis of existing data as those of the Cassini experiment [33]. In the future, it can also be reached by high 
accuracy Eddington tests (see for example the LATOR project [53]) or global mapping of deflection over the sky (see 
for example the GAIA project [54]). The former test would lead to an accuracy on the measurement of 7 sufficient to 
see the effect if its order of magnitude corresponds to that of the Pioneer anomaly. The latter test would allow one 
to reconstruct the dependence of 7 versus the impact parameter p and then, the dependence of <&p versus r. 

To summarize, the previous discussion of the effects associated with metric extensions of GR comforts the possibility 
to modify GR at large scales, while preserving the compatibility with existing classical gravity tests and allowing to 
include the Pioneer probes test. It also indicates that metric extensions could lead to new and observable effects at 
distance scales typical of the solar system. This provides further motivation for performing new experiments in the 
solar system in order to improve our knowledge of gravity laws, in particular of their scale dependence. 
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